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PREFACE 


The U. S. Naval Engineering Experiment Station, 
Annapolis, Maryland has recently run tests to determine 
the effects of thermal shock on steam piping for naval 
use. After several preliminary tests, data were taken 
for a section of pipe heated by the flow of saturated 3team 
through the pipe. These data consisted of temperature 
readings at the inside and outside walls of the pipe 
at time intervals of five seconds. The calculations 
of the stresses set up by thermal shock require a know¬ 
ledge of the temperature distribution in the wall at all 
radii. 

With the hope of assisting the Engineering Experiment 
Station in the evaluation of these test data, the author 
chose the task of finding the means of determining temper¬ 
ature distribution in a pipe wall at unsteady state con¬ 
ditions • 

The author wishes to express his deep appreciation 
to those who have aided in this investigation. He is 
especially grateful to Associate Professor George H. Lee, 

IT. S. Naval Postgraduate School, without whose assistance 
and continued guidance this investigation would not have 
been carried out and for his careful review of the complet¬ 
ed manuscript. 

The author is also grateful to the Engineering 
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Experiment Station and its Doctor William C. Stewart for 
suggesting the investigation and cooperating to the fullest 
in furnishing data. He wishes further to express apprecia¬ 
tion for the counsel of Doctor Warren M. Rohsenow, Mass¬ 
achusetts Institute of Technology, Cambridge, Massachusetts. 
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CHAPTER I 


INTRODUCTION 

A metal pipe initially of uniform temperature, is sub¬ 
jected to thermal shock by the sudden passage of a heating 
medium through the pipe. Y/hat the temperature distribution 
is in the pipe wall before steady state conditions are at¬ 
tained is the subject of this investigation. The author 
first became aware of this problem in August 1947 during 
an informal discussion with Doctor W. C. Stewart of the EES. 

Thermal shock in the operation of a Navy steam propul¬ 
sion plant is an abnormal condition. The possibility of 
thermal shock conditions through a casualty is however always 
present. High temperature shock can be the result of failure 
to properly "warm up" a main steam line before admitting high 
pressure saturated steam from the boiler. Cold shock can be 
the result of "carry-over" of a water slug from the boiler into 
the lines which have been carrying superheated steam. The 
problem was assigned the ESS for investigation. 

By August 1947 the EES had obtained a sizeable amount of 
data from cyclic thermal shock tests. These tests were run on 
various sections of pipe to determine the physical effects of 
cyclic thermal shocking. Heating had been accomplished by 
using exhaust gases, superheated steam, saturated steam and 
induction coils. Saturated steam produced the greatest temp¬ 
erature gradient through the pipe wall and therefore the 
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greatest shock. Cooling had been done by water or com¬ 
pressed air. 

Most of the tests indicated temperature readings only 
for thermocouples located at the inside and outside walls. 

A few test specimens were fitted with three thermocouples, 
the third thermocouple being located at the mean radius of 
the pipe wall. The reading of the thermocouples were taken 
at time intervals of five seconds. 

Calculations of thermal shock stresses had been made 
by the EES for representative tests. The author disagrees 
with calculations made by the EES for the determination of 
temperature distributions in the pipe wall in the unsteady 
state condition. Their means of solving for the temperature 
gradient is a graphical solution accomplished by a step-by- 
step method of averages for layers taken such that ig/r^ = 
r_/r =? a constant. The method fails to fit the data when 

X 8. 

the outside wall begins to increase in temperature. In 
the particular problem discussed in this paper that con¬ 
dition arises at time equal to six seconds. To allow for 
this inconsistency the EES progresses by adjusting the data 
to agree with the boundary conditions. 

Literature dealing with the analytic solution of 
temperature distribution in a hollow cylinder is limited 
to highly idealized problems. Peculiarly, all published 
solutions have one common boundary condition, namely, 
an instantaneous surface heat source which maintains the 
surface at a constant temperature. This condition is not 
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present in the EES problem nor in any other actual problem. 

The combination of three conditions is the cause for 
the complexity of this particular problem. These three 
conditions are: 1} steel pipe has a high thermal conduc¬ 
tivity; 2) the heat source in the EES problem produces a 
variable surface temperature and; 3) the heat source pro¬ 
duces a variable rate of change of the surface temperature. 

The analytic solutions published in the literature are of 
little value as related to the EES problem. 

Graphical solutions in the literature are developed 
for the elementary boundary conditions only. The theory 
of these, however, can be applied to the ESS problem. The 
author carries out this application for the case of heat- 
ing only and presents results which he feels are more accurate 
than those obtained to date by the EES. 

Attempts to find an analytic solution have proven un¬ 
successful. In the light of the authors investigations, 
there does not appear to be an analytic solution to this 
problem at present. 

It has become quite obvious, as a result of this 
investigation, that more detailed data should be taken in 
future tests. The installation of several thermocouples 
at proper radial points in the pipe walls would be desir¬ 
able. The maximum temperature rate and consequently the 
maximum thermal shock in the test occurs within five seconds 
after heating begins. Automatic recording devices attached 
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to the thermocouples may he the answer to better data in 
the critical time interval. 

The numerical results in this paper are only as 
reliable as the data taken. All values of temperature 
prior to a time of five seconds have been obtained by 
extrapolation. A plot of the data taken when three 
thermo.couples were installed in the test indicates a 
lag in the temperature readings,* This may be an error 
introduced by the use of thermocouples. The lag on the 
other hand may be only the result of personnel errors in 
taking the temperature readings. It is necessary there¬ 
fore to consider the value of this paper primarily in 
terms of the theory presented and its application to a 
specific problem. 

* C.f. Appendix V. 
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CHAPTER II 


ORIGIN OP THE DATA 

Tests were run on a length of 25-20 Cr-Ni, 3/4 inch 
IPS, schedule "double extra strong" pipe using saturated 
steam, superheated steam and exhaust gas as the heating 
medium. The resultant curves of inside and outside 
temperatures vs. time are shown in figures 1, 2 and 
3 respectively. 

One thermocouple was peened to the outside wall at 
mid length of the tube. Another thermocouple was brazed 
to the bottom of a £ inch diameter hole drilled to within 
0.01 inch of the inside wall, 1/8 inch distance from the 
peened in thermocouple. In the test with exhaust gas, 
however, the latter- thermocouple was pressed against the 
bottom of the hole by means of an 18-8 Cr-Ni steel screw 
insulated from the thermocouple. The thermocouple wire 
used was No. 29 B & S gage iron-constantan. The inlet 
and outlet pipes carrying the heating medium were drawn 
up tightly against the counterbored ends of the test 
sample by means of bolts secured to flanges welded to the 
pipes. Insulating gaskets were placed in the counter- 
bored ends of the test sample to minimize heat transfer 
due to conduction. 

The specially bladed turbo supercharger wheel at the 
EES was employed to supply the combustion gas using gas 
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temperatures ranging from 1040P to 1360F. Other runs 
were made by passing superheated steam at 500 psi pressure 
and 1050-1070F temperature through the sample. A third 
set of test runs was made by passing saturated steam at 
500 psi pressure and 470F temperature through the test 
pipe. 

In addition to these runs a series of tests was 
made on a similar sample of 18-8 Cr-Ni steel pipe, 3 1/4 
inch O.D. and 13/4 inch I.D. with thermocouples located 
.041 inches and .403 inches from the inside wall and 
a third peened to the outside wall. Thermal shock 
was created by means of induction heating. A plot of 
the temperature vs. time curves for all three thermo¬ 
couples is shown in figure 4. 

Runs were also made on an 18-8 Cr-Ni pipe, 31/2 
inches O.D. and 2 inches I.D. with thermocouples located 
.01 inches and .306 inches from the inner wall and a third 
peened to the outer wall. Thermal shock was created 
by means of exhaust gas at temperatures ranging from 
1260-1450F. A plot of the temperature vs. time curves 
for all three thermocouples is shown in figure 5. 

The results of these tests clearly indicated that 
saturated steam produced the greatest rate of temperature 
rise. Thermal shock tests were then planned using sat¬ 
urated steam as the heating medium. A section of two 
inch 25-20 Cr-Ni pipe 1.189 inches outside radius. 
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,749 inches inside radius was set up for test in the same 
manner as the initial tests. Thermocouples were located 
as for the previous test runs with saturated steam. A 
representative plot of the resultant temperature curve is 
shown in figure 6. 
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CHAPTER III 


ANALYTIC SOLUTION 

To calculate the thermal stresses in a finite hollow 
cylinder the following formulae given by Timoshenko (16)*' 
are used: 



It Is. not the author's intention to show the 
derivation of the above formulas since they follow a 
straight forward method and are presented clearly by 
Timoshenko'. The only purpose in introducing them is 
merely to show their relation to the temperature dis¬ 
tribution in a thick cylinder. 

To solve the integral Jq rdr it is necessary to 
obtain the temperature distribution through the pipe 
wall. It is therefore necessary to determine the 
analytic solution or to use a graphical method of 
solution. The Schmidt graphical method for a slab 

'^Numbers in brackets refer to bibliography. 
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( 14 )( 13 )( 11 ) can be adapted to a cylinder wall but large 
errors are introduced during the early stages of the 
solution. 

The maximum thermal stresses will occur in the early 
stages of heat transfer. It is therefore readily seen that 
an analytic solution is more desirable than the Schmidt 
graphical solution. 

1. The general solution. 

Assume the initial uniform temperature of the body as 
zero, the boundary conditions of the problem are: 


at 

o 

ii 

-p 
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r = r 

a 

• 
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t =0 

and 

r > r 

a 

• 

9 
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The differential equation for heat conduction in a 
cylinder (Fourier-Poisson) when radial changes only are 
observed is: 


(3) 


a q 

^ t 


- o(. 



2 + 


1 ^ e 

r <5 r 


A general solution is obtained by assuming Q=R T where R 
is a function of r only and T is a function of t only. 
Differentiating and substituting: 
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(4) 


2 


1 

oCT 


l T 
3 t 





o 

where^ is any arbitrary constant. 

A general solution therefore since T = e^ c l 

and R = A^J 0 (/<r) + B 1 N 0 ( y «r) 

is: 


(5) ©= ^ (AJ" 0 (^r)+BN 0 (^ r} e -0 ^ 

l 

where A and B are constants of integration. 

A particular solution occurs at steady state when 
<)©/,} t = 0 the solution having the form: 

(6) © = C log r + D 

where C and D are constants of integration. 

Under the conditions where the pipe wall eventually reaches 
a steady state condition the general solution is: 


(7) 

2. Solutions to similar problems. 

Barker (1) and Lees (9) give the steady state 
temperature distribution which is of interest but of little 
value• 

Barker’s form: 


© = C log r + D - 



J 0 (r) + BN 0 (/fr) e 


- 2 t 
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( 8 ) © 


9 a lQ S r/r b + © b log r a /r 
log r a /r b 


Churchill (3) and Williamson and Adams (17) present 
papers which give the analytic solution for the transient 
temperature distribution in a solid cylinder. The initial 
temperature is uniform. The surface temperature rises 
instantaneously to a given temperature and is maintained 
at that value. 

Williamson and Adams form: 

(9) 9=6, 1 - g \ -iqfeaAft-L "-vClA/r-bf* 

j A J l^n) 

— yn-l _ 

where yU is the nth root of J o (/0 = 0. 

Carslaw and Jaeger (2)* give a solution for the case 

of the infinite hollow cylinder. 

Boundary conditions: 

at t = 0 and r-r ; © = 0 

a 

at t > 0 and r=r j ©= ©. 

& JL 

Solution: 


(10) © = © 1+ ££i L-qC jj?t J 0 (^r)Y 0 (Mr a )- Y 0 (^r)J 0 ^r a ) d/A 

+ 4^ r o) ^ 


^C.f. Appendix I 
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Dahl (4) gives a solution for the case of the finite 
hollow cylinder. The initial temperature is uniform. The 
inner surface temperature is suddenly raised to a value 
and maintained at that value. The outer wall temperature 
is maintained at the initial temperature. 

Boundary conditions: 
at t=0 and r=r a j 0=©x 
at t - 0 and r > r a j 0 = 0 
at t = t and r = r a ; 0 = ©x 
at t = t and r=r b ; 0=0 
Solution: 


n=°o 


( 11 ) G = 


log r^/rj 



n=l 


- J o 


r o^*n r ^“ „ v 

,I o ( /'n r b ) 


»o(**n*> 


e 




(lla) 



log r b /r a 


Jn(Mnrb) 

N o^n r b) 


N l</*n r a>~ J l^n r a) 


(lib) A’o = 


2 ©i/log rt/r a 



J °—— Jx(^ n r b ) 
N o^n^) 




2 0 X / log rh/r a _ 

J lOVa 

K'o'/'a^b) 
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ts to 



Where jul has the values of the positive roots oft 


(11c) 


J o <A r a > 


Jp Wn r b J 

No ^n r b ) 


VAiV 



Jahnke und Emde (6) present excellent tables for the 
evaluation of both jx and the first two orders of N. 

Carslaw and Jaeger (2) also present the case of a 
finite hollow cylinder. The inside surface temperature 
is raised instantaneously and maintained constant. The 
outer wall is adjacent to an infinite fluid which maintains 
it3 temperature constant at a point distant from the 
outer wall. 

Boundary conditions: 

at t =0 and r^r^r^j 9=0 

at t > 0 and r = r fi j 9 = 9-^ 

at t>0 and r = r^ ; 59/^r + f/k (9)=0 

where f is the unit surface conductance. 

Solution: 


(12) 9=9 X K 



. T h J o C^nPb) ~ /*n J 1 CMn*»b) 1 
1 


[ J o(Mn*> Y o^n*a> 

(k + /n^ , 3 o^n r a^ 


K(A r b) -A J i(Abb) 
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1 + r, 

_b 

1 +r b 


h log r^r 
h log r b/ r a 


where jJt are the positive roots of: 


(12a) J 0 (/^r a ) 


llY o ( y wr b ) ‘ y MY l ( > (r b ) J “ 

'*ol/* T a) [ h;r o ( / Ar b ) -^ J i^ r b)] 


= 0 


where h is the ratio f/k. 

3« Further investigations. 

Dahl’s equation appeared to be applicable at first. 

Upon further investigation it was found that the actual bound¬ 
ary conditions were not in agreement with Dahl’s boundary 
conditions. An attempt to apply temperature differences in a 
step-by-step manner failed at the second step since Dahl’s 
equation no longer applied without an initial uniform temp¬ 
erature distribution. Further consideration of Dahl's 
equation was therefore abandoned. 

For any one radius in the EES problem © = 0(t) is of 
the form of Froelish’s equation (13) for the saturation 
curve. 


(13) 


0 = 


x t 
y + t 


where x and y are arbitrary constants to be determined from 
the shape of the curve. 
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At t = °° and r = r ; 9 = © 

a a 

at t=o° and r-r^; 9 = 9 fc 

at t = °° and r a <r<r b ; ©^Barker’s steady state distribu¬ 
tion* If the general solution were: 

c»o 

(14) 9 = C + ^ g.J o (y^r) + BN Q (y«r)] e -*>« 2 t 

i 

then at t~°° the second term disappears. Let C equal the 
steady state condition and the final boundary condition 
is satisfied. Although both inside and outside wall temp¬ 
eratures are functions of time, they also satisfy this 
boundary condition. They become constant at infinite time. 
At time equal to zero and r= r^ the Barker term becomes 
zero (initial temperature) and the constant B of the Bessel 
term can be evaluated. Also at time equal to zero and 
r =r a >y< must satisfy the roots of the Bessel term in the 
above equation. It now remains only to evaluate the con¬ 
stant A. This unfortunately is impossible. At zero time, 
zero temperature, r-r, there exists only the value of zero 
for A and the solution thus fails. 

This practical problem fails to produce the very con¬ 
venient boundary condition of an infinite surface heat 
source. The EES is working with iron pipe of very high 
thermal conductivity. The problem primarily is centered 
around the initial two or three second period of time. 
Boundary conditions used by Dahl (4) or Carslaw and Jaeger 
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(2) cannot be assumed. The problem cannot be stretched as 
in the case of Perry and Berggren (12) who carefully choose 
the cork lagging on a metal pipe and then work only with 
the cork cylinder. 

Another particular solution is to let ^©/^t=C 1 . 

The solution of the differential equation is then:* 

g 

(15) —+B 1 In r t + Cj_t 


Where and are constants of integration; 

Another general solution therefore might be: 


(16) ©-■ 


C^r* 


■tB^ln r-t-D^-f-C^t 


-W Wr , + BNoH 




This equation can be used only for a limited time since the 
first term Increases with time. The entire solution breaks 
down on a substitution of initial conditions. At all values 
of r, time equal to zero, the first term and the temperature 
difference (© a - © b ) are both zero. Again A equal to zero 
Is the only choice for evaluating A.. Unless © is considered 
as © a ~ ©^, the solution is merely for a condition in which 
all temperatures are rising at the same rate at all radii. 

For a perfectly insulated outer wall the problem is 
only slightly modified.. The outer and inner walls rise to 


* C*f. Appendix II 
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the same value at infinite time. The Bessel term of the 
general solution for the differential equation is alone 
not enough since it becomes zero at infinite time. Assume 
a general form: 


(17) r,t) - 



+■ BN, 




e 


-oC/^t 


I 

where $(r,t)-Cg when t=°° 
and j6(r,t) = 0 when t = 0. 

Boundary conditions again permit evaluation of B and />> 
hut again all terms go to zero when we attempt to evaluate 
the constant A. 

Assume the solution of the differential equation to 
be ©=R+T. Differentiation and substitution results in: 


(0-8) 


1 . ( 2 2 R . 1 >R \ _ ,.2 

“it 


On inspection it is seen: 


(19) 


l^-scu 2 = Je 
^t ^t 


thus the result is the same as the solution of ^©/^t-C^* 

As a result of the above investigation there does not 
appear, applying the above approach, to be an analytic 
solution to the EES problem. Through necessity, rather than 


17 



through preference, an attempt will be made to find a 
graphical solution* 
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CHAPTER 17 


GRAPHICAL SOLUTION 
1. Available methods. 

Schmidt’s graphical method (14) is basically for a 
slab. It can be adapted to the cylinder, but as previously 
stated, it is very much in error during the early time in¬ 
tervals. Williamson and Adams (17) point out, and Bishenden 
and Saunders (5) confirm the fact, that Lord Kelvin’s (10) 
probability integral more nearly approximates the initial 
thermal wave penetration. This improvement over the Schmidt 
method was further advanced by Perry and Berggren (12). 

They determine a graphical solution for the heat flow 
differential equation for cylinders using small finite 
differences.* This graphical solution is developed only 
for the problem of a hollow cylinder with an instantaneous 
surface heat source. The inside wall temperature is then 
maintained constant. Conduction to an infinite fluid is 
assumed to take place at the outer wall. The results 
are compared with values obtained by using Carslaw’s 
and Jaeger’s formula for the same condition.** 

The equation given by Perry and Berggren is for the 
temperature rise during a time interval at any given 

* C.f. Appendix III 

**C.f. Chapter III, equation (12). 
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radius (12) (7). It is merely a weighted mean value of 
the temperature of the two adjacent radii before the time 
Interval, namely: 


(20) 9 


t + At,T 




+ (!+-££ 


2r 


-)®t. 


r + a r 


For very small values of at or as dr-?-0 the equation 
agrees with Schmidt’s equation for a slab (14): 


( 21 ) 


+ -a t 




,r- nr 


+ 9 


t ,r + a 



The initial curve following the Lord Kelvin probability 
integral or Dahl’s equation* could be used over that first 
small time increment. Assume an instantaneous rise of 
temperature at the inner wall for the short time increment. 

The outer wall will remain at the initial uniform temper¬ 
ature. Where dt is very small Dahl’s equation has little, 
if any, advantage over Lord Kelvin’s integral. 

2. A graphical solution of the EES problem. 

It is desirable to have small increments of dr. 

2 ** 

Since At is proportional to dr , the error of the graphical 

* C.f. Chapter III, equation (11) 

** C.f. Appendix III. 
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solution can be reduced to a minimum by a proper choice 
of Ar. By progressing through the graphical analysis 
using a tedious mathematical solution for each point, an 
accurate curve can be drawn. In the EES problem the outer 
surface temperature increases with the inside temperature. 
It is desirable to prevent a dip in the curve below the 
outer wall temperature. Any temperature which graphically 
is determined to be below the outer surface temperature 
will automatically be given the same temperature as the 
outer surface. This condition will occur during early 
time intervals near the outer layer. An attempt to fair 
in a curve would introduce possible errors on the unsafe 
side of the actual temperature distribution curve. It 
is considered better therefore to introduce a slight 
safety factor by using the suggested procedure mentioned 
above. A similar* procedure is used by Sherwood and Reed 
(15) for the slab. 

Jakob (7) and Perry and Berggren (12) give proof 
that a scale of Ar plotted as A A, where ,A= In r/r Q , 
can be used in a plot against temperature to produce 
a direct graphical method of solution.* The temperature 
rise in a cylindrical segment for a time interval At 
can be obtained merely by drawing a straight line on 

the graph between -A A 6111(1 ®t,l+-AA * The 

intersection ©t + A t, 7\. is tlie desired temperature. 

* C.f. Appendix IV 
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A careful evaluation of the problem graphically 
yields the results as shown in figure 7 and figure 8. 
It appears, therefore, that the graphical solution 
suggested by Perry and Esrggren as applied to the EPS 
problem is the only solution if it is desired to main¬ 
tain the complex boundary conditions* It is also a 
solution which can be calculated with comparative ease 
by a direct graphical method. 
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APPENDIX I 


Carslaw and Jaeger (2) present the following solutions 
in their paper "Some two-dimensional problems in conduction 
of heat with circular symetry”. 

Cylinder of radius a, initial uniform temperature 0, 
surface temperature suddenly brought to 6-j. and maintained 
constant. 

3. Infinite hollow cylinder of radius a, initial uniform 
temperature 0, inner surface temperature suddenly brought 
to ©2^ and maintained constant. 

5. Solid cylinder of radius a» initial uniform temperature 

©, and radiation at the surface takes place into a medium 
» 

at temperature 0. 

6. Infinite hollow cylinder of radius a, initial uniform 
temperature ©^ and radiation takes place into a medium at 
temperature 0. 

7. Hollow cylinder of initial uniform temperature 0. 

a) Outer surface maintained at 0 while inner surface 
is suddenly raised and maintained at ©}_. 

b) Inner surface maintained at 0 while outer surface 
is suddenly raised and maintained at ©^. 

c) Inner surface raised suddenly to ©„ and maintained 
while outer surface is suddenly raised to ©^ and 
maintained, 

8. Hollow cylinder of initial uniform temperature 0, 

i) Outer surface suddenly raised to 9v> and maintained 
constant while inner surface radiates to a medium 
of temperature © a . 

* Numbers of problems correspond to those used by Carslaw 
and Jaeger. 


25 



ii) Inner surface suddenly raised to 9 and maintained 
constant while outer surface radiates to a medium 
at temperature G^, 

iii) Inner surface radiates to a medium at temperature 
9 a while outer surface radiates to a medium at 
temperature G^. 

9, Solid cylinder with an instantaneous heat source at 
the axis, 

a) Outer surface temperature maintained at 0, 

b) Outer surface perfectly insulated, 

c) Outer surface radiates to a medium at temperature 0. 

10, Solid cylinder with an instantaneous surface heat 
source at any radius r, 

a) Outer surface temperature maintained at 0. 

b) Outer surface perfectly insulated, 

c) Outer surface radiates to a medium at temperature 0. 

11, Solid cylinder with an initial temperature a 0(r), 

i) Outer surface temperature maintained at 0. 

ii) Outer surface perfectly insulated, 

iii) Outer surface radiates to a medium at temperature 0. 

12, Infinite hollow cylinder with an instantaneous surface 
heat source at any radius r, 

a) Inner surface temperature maintained at 0, 

b) Inner surface perfectly insulated, 

c) Inner surface radiates to a medium at temper¬ 
ature 0. 

13* Infinite hollow cylinder with an instantaneous surface 
heat source at any radius r and initial temperature a 0(r). 
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i) Inner surface temperature maintained at 0. 

ii) Inner surface perfectly insulated, 

iii) Inner surface radiates to a medium at temper¬ 
ature 0, 

14* Hollow cylinder with an instantaneous surface heat 
source at any radius r, 

a) Inside and outside walls maintained at temper¬ 
ature 0* 

i) Outer wall maintained at temperature 0 while the 
inner wall radiates to a medium at temperature 0, 

ii) Inner wall maintained at temperature 0 while the 
outer wall radiates to a medium at temperature 0. 

iii) Inner and outer walls both radiate to a medium at 
temperature 0« 

15. Hollow cylinder with an instantaneous heat source at 

any surface of radius r and an initial temperature a 0(r). 

i) Inside and outside walls maintained at temper¬ 
ature 0, 

ii) Outer wall maintained at temperature 0 while the 
inner wall radiates to a medium at temperature 0, 

iii) Inner wall maintained at temperature 0 while the 
outer wall radiates to a medium at temperature 0, 

iv) Inner and outer walls both radiate to a medium 
at temperature 0, 
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APPENDIX II 


A particular solution to the differential equation 
for heat flow in a cylinder can be developed for ^©/H = a 
constant. The differential equation becomes: 

d 2 ©/dr 2 + l/r . d©/dr =C 1 
Let d©/dr=y, then 

dy/dr + l/r . y -C]_ 

Let v=y/r, then 

dy/dr = v r dv/dr 

or 

2v + r dv/dr 
dr 

- .f 

r 

In r + Cg + 1/2 In (C^- 2v) = 0 
(r 2 ) (C -2v) = B 

Since v^y/r 

C^r 2 - 2 yr = B 

Since y = d©/dr 
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C-jT 2 - 2r d©/dr = B 
Jc-^r dr -J 2 d© = B J "■§ ? ,.. 


Cjr* - 4© = B» In r 4-D 


© - £ C x r^ + B* f In rtD' 

A particular solution to the partial differential heat 
transfer equation: 


( 3 ) rf f) 2 * t 1 ) £ 

^ r 2 r <) r 


J©_ 

) t 


is 


(15) 


© = £ C-jr 2 + B x In r + d 2 t Cjt 


when 


^©/^t - a constant Cj 
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APPENDIX III 


A graphical interpretation of the differential 
equation for heat transfer in cylinders is presented by 
Ferry and Berggren (12)* The differential equation 
when axial temperature changes are zero is: 

(3) -2i_ = oc * _i_ 1 ®-\ 

it (^r 2 r ir J 

In terms of small finite differences the above equation 
becomes: 



rt-flA. A A+4A 


Figure 9 

Basis for a graphical solution 
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From figure 9 the following equations can be readily- 
determined : 


A 9 _ ^t + a t ,r *~ ®t ,r 
At At 


A 9 9 -9^ 

__ t,r + a r t, r-Ar 

Ar 2 A r 

9 _ g 9 9 

A 2 e _t i rj_A_r -£jT -t^Ar 

_ -_ A£ _ A£ _ 

Ar 2 Ar 


^t,r+ at " 2 ^t,r + ^t,r~ nr 
(Ar) 2 

Upon substituting these values into the equation of 
finite differences: 


JL_ 

oC 



a t,r~ Q t,r 


At 


Q t,rt a r-^^t.r't e t,r- 
( Ar) 2 


•A r 


1 Q t ,r + a r*~ e t ,r- a r 
r 2 Ar 


Let At = ( A r) 2 /2cA , and multiplying by ( A r) 2 : 
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2(©. ^ -e 

t+4t,r t,3 


dr 

2r 


^t ( r+flr ^t # r-i5r) 


(20) e t +at ,r^ 



-) e 


t ,r+ hr 


+ ( 1 - 


£2L.)e. 

2r t,r- Ar 


The resultant equation is one which is a weighted mean 
favoring the side of increased radius and the side of 
greater heat capacity. 
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APPENDIX IV 


A plot of temperature vs. In r/r can be used for a 

6L 

direct graphical method of solution. 



Figure 10 

Scale for the direct graphical method 
From the geometry of the above figure it is apparent 

that: 


®t + a t ,r 


AA 2- «t 


■ r- ar + Vtt a r 

^ A^ + AA 2 


But from 


(20 > e t + . t,r* * [ (1 + g> e t,r + . r + U- §)e t , r . 
AAt = i+*Z AA~ -i- 


hr 






AA 2 

AXi 



r + 


AT 

2 


but r - Ar/2 and r+nr/2 are the mid radii of the two 
adjacent Ar layers of r. The scale division therefore 
must be inversely proportional to their corresponding 
radii. 

This is satisfied by choosing a scale variable 

A = In r/r Q 

so that dA = dr/r or AA~nr/r 

The mean radius of AA^ is r - Ar/2 and similarly for 
A A 2 is r + nr/2, 

Then: 

a „ nr _ _Ar_ 

r ■. flr t* j- Ar 

2 2 


and 

AA 2 _ r- ± 

AX-, r ., Ar„ 

2 

These value coincide with the values required by figure 10, 
This scale for plotting radius therefore provides a means 

of making a direct graphical solution. 
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APPENDIX V 


A plot of steady state temperature distribution radial¬ 
ly in a hollow cylinder is a straight line when plotted on 
coordinates of temperature vs. In r/r a . The data as shown 
in figure 5 represents the readings of three thermocouples 
in a pipe wall. This data was replotted in figure 11 as 
temperature vs. In r/r a . 

Two peculiarities are apparent upon examination of the 
resultant curves. Most apparent is the change of curvature 
from concave up to concave down after ninety seconds. The 
second peculiarity at ninety seconds is a change from increas¬ 
ing slopes with time to decreasing slopes with time. The 
latter is the effect of the wave front as the heat penetrates 
the pipe wall. 

The change from concave up to concave down indicates 
at once a descrepancy in the data. Since steady state is 
the ultimate condition which can exist under the conditions 
of the test, the concave down curves are an impossibility. 

The author is of the opinion that the plot reveals a 
lag in temperature readings. This may be an error introduced 
by either the thermocouples, the indicating instrument, or 
both. The lag, on the other hand, may be only the result of 
personnel errors in taking the temperature readings. The 
existence of this lag must, however, be acknowledged. 
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APPENDIX VI 


The tedious graphical solution represented by figure 
7 immediately caused the author to investigate the possibil¬ 
ity of a more rapid graphical solution. A graphical solu¬ 
tion, represented by figure 12, was then made for purposes 
of comparison. 

The first sixteen steps of both solutions are ident¬ 
ical. Time increments in figure 12 are then rapidly in¬ 
creased to speed up construction. The time to produce 
figure 12 was approximately one half the time required for 
figure 7. A comparison of resultant values of temperatures 
for r]_, rg and r§ from both figures was then made. 

A maximum deviation of 6 degrees representing a 4$ 
error occured for r^, at 1.512 seconds. This error rapidly 
decreases and disappears after 6.048 seconds. For rg 
the maximum deviation was 4 degrees for a 2.8$ error at 2.1168 
seconds. The trend of a decreasing error at an increasing 
time for that maximum error with an increase in radius was 
noted for rg. At rg the maximum deviation was 3 degrees 
for an error of 2.7$. The time was 3.628 seconds. The 
deviation tends to disappear at r 2 and r 3 as time increases 
this is a trend characteristic of any graphical solution of 
this type (15). 

The maximum error introduced by the rapid graphical 
solution is small and tends to decrease with time. Deter¬ 
mination of the temperature distribution in a pipe wall 
can be made by the rapid graphical method in two to three 
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hours. These factors should be of value to the EES who 
may wish to make calculations for a series of thermal shock 
tests for comparative purposes. 
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